Off-diagonal disorder in the Anderson model of localization 
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We examine the localization properties of the Anderson Hamiltonian with additional off-diagonal 
disorder using the transfer-matrix method and finite-size scaling. We compute the localization 
lengths and study the metal-insulator transition (MIT) as a function of diagonal disorder, as well as 
its energy dependence. Furthermore we investigate the different influence of odd and even system 
sizes on the localization properties in quasi one-dimensional systems. Applying the finite-size scaling 
approach in conjunction with a nonlinear fitting procedure yields the critical parameters of the MIT. 
In three dimensions, we find that the resulting critical exponent of the localization length agrees 
with the exponent for the Anderson model with pure diagonal disorder. 



Introduction. The electron localization properties of disordered solids are of high interest both experimentally and 
theoretically. A simple approach to this problem is given by the Anderson model of localization. In this model one 
considers a single electron on a lattice with N sites described by the Hamiltonian 

N N 
i^j i 

where \i) denotes the Wannier states at site i. Disorder is usually introduced by varying the onsite potential energies 
Ci randomly in the interval [—W/2,W/2\. In this work we report on the effects of additional off-diagonal disorder 
i.e., with random hopping elements tjj £ [c — w/2, c + w/2] between nearest neighbor sites. Thus c represents 
the center and w the width of the uniform off-diagonal disorder distribution. The energy scale is set by keeping w = 1 
fixed. This kind of disorder is similar to the random flux model ^|,|5| , but here the modulus of the ty is random while 
the phase is constant. The Hamiltonian ([l]) has a chiral symmetry on a bipartite lattice, i.e., under the operation 
ip n — > (— \) n tj) n , the Hamiltonian H changes it sign and as a consequence, if e„ is an eigenvalue of H , then so is — e„. 

We calculate the localization length A of the electronic wave function by the transfer-matrix method (TMM), which 
is based on an iterative reformulation || of Eq. (|l|) as 

</Wl \ =Tn (^ \ = ( *n+l( S1 - H i) -Klltn \ ( V>» \ (2) 

where E is the energy, Hj^ is the Hamiltonian in the nth slice of the TMM bar, 1 and are unit and zero matrix, 
respectively (7), and the diagonal matrix t„ represents the hopping elements connecting the nth with the (n — l)th 
slice. The cross section of the bar is 1, M, and M 2 for spatial dimension d = 1,2, and 3, respectively, and we always 
assume periodic boundary conditions for the TMM. Starting with two initial states (ipi,i>o) one iterates along the 
quasi ID bar until the desired accuracy is obtained. The localization length = l/7 m i n is computed from the smallest 
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Lyapunov exponent 7 m i n of the eigenvalues exp(±7i) of lim„^ 00 (r T ^r„) 1 / 2 ™ with r„ = T„T„_i • • • TiT\. Assuming the 
validity of one-parameter scaling close to the MIT, we expect the reduced localization lengths Ad(M) /M to scale onto 
a scaling curve \d{M)/M = fd{£,/M) with scaling parameter £. 

Localization lengths at the ID Dyson singularity. In ID it has been shown that the density of states (DOS) of 
the disordered model diverges at E = |||| , if W = on a bipartite lattice || . This singularity is intimately related 
to a divergence of the localization length ||[l0| . Assuming that the localization properties of the wave function can 
be described as usual by an exponential decay, one finds analytically 



(3) 



with vp the Fermi velocity and D an energy parameter depending on the strength of the disorder. In Fig. [I] we show 
that this energy dependence is convincingly reproduced in numerical data of the model: Ai diverges at E = as in 
Eq. (||) for W — > 0. Note that the smallest values for the localization length are obtained for c = 0.25. We remark 
that the validity of the exponential form of the wave function has been questioned in previous studies |3|,|l2]| presenting 
evidence for power-law localization usually related to critical states. 

Odd/even effects in quasi ID systems. In 2D, the DOS also has a sharp peak at the band center |l]-|j| and it 

was shown by the TMM that the localization length again diverges HQ. A finite -size-scaling (FSS) analysis of the 
TMM data together with studies of the participation numbers and multifractal properties of the eigenstates revealed 
|p3| , H that the states at E = show critical behavior like in the 3D Anderson model at the MIT [[l5) . 

Renewed interest in the study of quasi ID systems with off-diagonal disorder stems from the fact that the functional 
form of the divergence depends on whether M is odd or even as shown by a Fokker-Planck approach in Ref. |l0| . This 
question is intimately linked to the presence or absence of bipartiteness for the given lattice and boundary conditions. 
In Ref. numerical data in support of the analytical results had already been shown. Here we will broaden the 
investigation by studying how the localization length at E = is influenced by the odd/even effects. Eilmes et 
al. |p^ , p^| observed that the reduced localization length A2/M for the 2D Anderson model with random hopping is 
independent of the system sizes at E = for bipartite lattices up to M = 200. Motivated by the results of Refs. ||[l0) , 
we show the behavior of the localization lengths for non-bipartite systems with odd M as well as bipartite systems 
with even M in Figs. I and § We see that the values of the reduced localization length differ for odd and even M. 
Nevertheless, A2/M at E = remains constant and thus critical in the sense of Ref. |ll| for odd and even M up to 
M = 101 and 180, respectively. 

Critical exponents in 3D. For the 3D model we computed the critical parameters at the MIT when either E or 
W are varied across the transition, i.e., £ cx \E — E C \~ VE or £ oc \W — W c \~ Uw |l6| at fixed c. The left panel of Fig. || 
shows the reduced localization length A3/M for even M up to 14 at E = and c = 0. The transition from extended 
states for W < 4.05 to localized states for W > 4.05 is clearly visible as A3/M decreases (increases) with increasing 
M for the extended (localized) case. In the right graph the resulting scaling function /3^/M) is plotted. It was 
obtained by a non-linear fit taking non- linear and non- universal corrections to FSS fl7|l into account. Comparing the 
exponents ve = 1.61±0.07 and vyy = 1.54±0.03 of the transitions obtained from the fits at W = 0, c = and E = 0, 
c = 0, respectively, with recent results for the Anderson model with pure diagonal disorder we find good agreement. 
Conclusion. In this work we have studied the localization properties of the Anderson model of localization with 
off-diagonal disorder. We find non-localized states in ID and 2D only at the band center. In quasi ID we examined 
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odd/even system size effects for the localization length. We showed that the wave functions at E — in 2D remain 
critical up to the system sizes considered regardless of the odd/even effects. Our numerical results match reasonably 
well with the predictions |||l0|]. In the 3D case, we obtained the critical exponents of the MIT using TMM and 
FSS. Their values agree with recent results for the model with only diagonal disorder jl7],|l1|. Thus although the 
off-diagonal case exhibits some unusual features, its physics is nevertheless accurately described by the orthogonal 
universality class within the scaling theory of localization |(|[l9| . 
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FIG. 1. Localization length Ai as function of energy E with off-diagonal disorders c = (□, ■), 0.25 (o, •), and 0.5 (A, A) 
for diagonal disorder W = 0.01 (filled symbols) and 0.1 (open symbols) obtained by TMM with 1% accuracy. Solid (dashed) 
lines indicate data for E > (E < 0). 
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FIG. 2. Reduced localization lengths A2/M as a function of energy for odd system sizes (left panel, M = 3, 7, . . . , 35 from 
top to bottom) and even system sizes (right panel, M = 4, 6, 10, 14, 18, 20, 22, 26, 30, 34 from top to bottom) with random 
hopping for c = 0.25. Note the strong decrease in A2/M at E = for the bipartite even case. 
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FIG. 3. Reduced localization length X2/M vs. 1/M for E = and even M up to 180 with 1% accuracy correspond- 
ing to c = 0, 0.05, . . . , 0.25, . . . 0.6 as indicated by a, b, . . . , f, . . . , m The circles represent data for odd system sizes 
M — 11, 13, . . . , 89 and M — 101 with c = 0.25, i.e., the smallest X2/M values, and 0.5% accuracy. 
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FIG. 4. Left: Reduced localization length A3/M as function of diagonal disorder W at c = and E — obtained by TMM 
with 0.1% accuracy. Note the MIT at W c » 4.05. Inset: Enlarged region close to the MIT, the error bar is shown separately. 
Right: FSS curve for the TMM data from the left panel. 



